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Abstract: 

We study 't Hooft anomaly matching in lattice models with strong Yukawa or multi-fermion 
interactions. Strong non-gauge interactions among the mirror fermions in a vector like lattice 
gauge theory are introduced with the aim to obtain, in a strong-coupling symmetric phase, 
a long-distance unbroken gauge theory with chiral fermions in a complex representation. We 
show how to use exact lattice chirality to analyze the anomaly matching conditions on chiral 
symmetry current correlators at finite lattice spacing and volume. We perform a Monte 
Carlo study of the realization of anomaly matching in a toy two-dimensional model with an 
anomalous niirror-fermion content at strong mirror Yukawa coupling. We show that 't Hooft 
anomaly matching is satisfied, in most of the phase diagram, via the minimal solution in 
either the massless fermion or "Goldstone" mode, while in some cases there are extra massless 
vector like mirror fermions. The mirror spectrum at strong coupling is thus consistent with 
long-distance unitarity. We discuss the implications of our results for future studies of the 
most interesting case of the decoupling of anomaly- free mirror-fermion sectors. 
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1. Introduction, summary, and conclusions 

1.1 Motivation 

Understanding the dynamics of strongly-coupled chiral gauge theories is a long-standing prob- 
lem, where little has been learned since the 1980s [l]."*^ While not motivated by a pressing 
need to explain current data, this issue has its own theoretical interest and appeal. Strong 
chiral dynamics may gain physical applications in the future, should results from the LHC 
indicate that strong interactions play a role in physics near the TeV scale. 

This paper is on a topic even more theoretical than that of studying strong chiral gauge 
dynamics — the problem of constructing a gauge invariant lattice formulation of chiral gauge 
theories. Different approaches and the status of this subject are reviewed in [5, 6, 7]. 

1.2 The idea 

The idea that will be pursued here is, in its essence, a relatively old one, to the best of 
our knowledge due to Eichten and Preskill (EP) [8] (but see also the work of Smit [9]). EP 
proposed to use the strong-coupling symmetric phases [10, 11, 12, 13, 14], which occur in 
lattice theories with non-gauge — e.g., four-fermi, or Yukawa — interactions, to decouple the 
mirror fermions from a vectorlike gauge theory. The goal is to thus obtain, at long-distances, 
an unbroken gauge theory with massless chiral fermions in a complex representation. 

The strong-coupling symmetric phases are a lattice artifact, similar to high-temperature 
phases in statistical mechanics, where correlations range at most over a lattice spacing with 
no interesting dynamics at long distances.^ The idea of EP, appropriately rephrased to suit 
our current context, is to let strong non-gauge interactions act only on the mirror fermions 
and use the short-range correlations to decouple them from the infrared physics. One can 
describe the decoupling as the binding, due to the strong interactions, of all mirror chiral 
fermions into gauge-singlet or vectorlike composites with mass of order the lattice cutoff. 
In this picture of strong non-gauge mirror dynamics, due to asymptotic freedom, the gauge 
interactions play only the spectator role of weakly gauged global symmetries (for a recent 
study of the decoupling of scalars in a strong-coupling symmetric phase in a simple toy model 
with gauge fields, see [15]). Ideas closely related to that of EP have received attention in the 
past. In all cases studied, it was found that either the mirrors fail to decouple from the long 
distance physics or the purportedly light fermions also obtain mass [16, 17, 18, 19, 20, 21, 22]. 

Despite the failed early attempts, a proposal to revisit the EP idea in light of the rela- 
tively recent progress in exact lattice chiral symmetry was made not long ago [23] (a suggestion 
along similar, but not identical, lines had been made earlier in [24]). This renewed interest 

^With the exception of some progress in chiral gauge theories with f\f=l supersymmetry [2, 3] and recent 
advances in the study of semiclassical chiral dynamics [4] . 

^ Often, the question of the "triviality" of Yukawa interactions is raised when strong Yukawa interactions 
are mentioned. We note that in our context, triviality — interpreted as the absence of any long-distance physics 
in the strong-coupling symmetric phase — is to be advantageously exploited, since the goal of the Yukawa 
interactions is to keep all mirror fermions at the cutoff scale (see, e.g. [8]). 
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was motivated by the fact that in the pre-1997 formulations of vectorUke gauge theories on 
the lattice chiral symmetries were explicitly broken. The absence of an exact chiral sym- 
metry without doublers makes it impossible to split a lattice Dirac fermion field into chiral 
components — one to be henceforth called a "mirror" and the other "light." As a result, any 
strong non-gauge interaction inevitably couples to both the "mirror" and "light" components 
of the vectorlike fermion. The analysis of the strong dynamics — clearly not purely "mirror" — 
of the chiral symmetries and their realization then becomes involved and ambiguous. Most 
importantly, previous formulations lacked a manifest symmetry explaining why some chiral 
fermions should stay exactly massless while others could obtain mass; instead, chiral sym- 
metries distinguishing "light" from "mirror" modes and protecting the "light" modes were 
expected to somehow emerge at special values of the couplings. 

In contrast, the recently discovered exact lattice chiral symmetries without doublers, 
defined via the Ginsparg- Wilson (GW) relation [25] and the Neuberger-Dirac operator [26, 
27, 28, 29, 30], permit the splitting of a lattice vectorlike fermion into "light" and "mirror" 
components. Ward identities for anomalous [31, 32, 33, 34, 35] and anomaly-free global 
symmetries hold at finite lattice spacing and volume, ensuring that the chiral symmetries 
that protect the light fermions from acquiring mass are exact. Thus, EP-like ideas for the 
decoupling of the mirror fermions can now be elegantly formulated using exactly chiral lattice 
fermions [23]. This elegant realization of symmetries comes at a price, however, as one still 
has to show that the strong mirror dynamics decouples the mirror fermions. The study of the 
mirror theory dynamics is now complicated by the exponential-only^ locality of the Neuberger- 
Dirac operator — this is in contrast with local-fermion formulations (whose drawback is the 
already-mentioned lack of a manifest "light"- "mirror" split) where the mirror dynamics can 
be studied via a relatively straightforward strong-coupling expansion. 

The lack of an obvious controlled expansion implies that the mirror fermion dynamics 
in a formulation of the EP idea with GW fermions has to be investigated via Monte Carlo 
methods. Due to the high cost of simulations with exactly chiral fermions, such studies are in 
their infancy. On the other hand, since the issues considered are ones of principle, analyzing 
two-dimensional models is a sensible first step. Clearly, showing that the ideas work in 
two dimensions will not prove that four dimensional chiral gauge theories can be similarly 
formulated, but the results are likely to provide insight into the relevant mirror dynamics. 
Another appropriate simplification is to neglect the gauge field fluctuations, since gauge fields 
play a spectator role to the strong mirror dynamics. If the strong non-gauge mirror dynamics 
gives the mirror fermions mass of order the lattice cutoff (in a manner alluded to above), 
asymptotic freedom of the gauge interactions and the exact chiral symmetry of the light 
fermions lead us to expect that turning on dynamical gauge fields will not significantly affect 
the mirror dynamics or lift the massless modes. Including dynamical gauge fields would thus 
only make sense after the mirror decoupling in zero gauge background is demonstrated. 



^Exponential locality holds subject to the "admissibility condition" on the gauge field background [36, 37]. 
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1.3 Its status and a question 

The one existing Monte Carlo study, by Giedt and one of us (E.P.) [38], is of tlie two- 
dimensional massless Schwinger model with strong mirror interactions added, a la EP, in an 
attempt to decouple one of the chiral components of the Dirac fermion. On general grounds, 
one expects that it should not be possible to decouple the mirror and have an unbroken 
gauge symmetry in this model, since the resulting long-distance theory — the chiral Schwinger 
model — would be anomalous [39, 40]. Nevertheless, this potentially "pathological" setup is 
the easiest allowing to address questions such as the phase structure of the mirror interactions 
with GW fermions, because of its simplicity and associated low cost of simulation 

In [38] , the strong dynamics of the mirror sector of this model was studied via a Monte 
Carlo simulation, in a vanishing gauge background. It was found that the (would-be-gauged) 
chiral symmetry of the mirror theory is unbroken"^ and, thus, the strongly-coupled symmetric 
phase exists in this model, in a large region of the mirror-coupling space. More intriguingly. 
Green's functions of local operators probing the charged fermion spectrum, elementary or 
composite, showed no evidence for long-range correlation — and thus for massless charged 
particles. These results are clearly puzzling, as gauging the unbroken global chiral symmetry 
would give rise to an anomalous unbroken gauge theory. 

The purpose of this paper is to gain a more detailed understanding of the resolution 
of this puzzle. In particular, we would like to know what happens when one attempts to 
decouple mirrors in anomalous representations using non-gauge strong interactions. There 
appear to be two possibilities, in a Euclidean rotation-invariant theory: a.) taking the mirror 
interactions of GW fermions strong leads to an inconsistent theory in the continuum limit^ 
or b.) there are massless charged states — Goldstone bosons or fermions — canceling the light 
fermion anomaly, which were somehow missed in [38]. 

We first addressed this question in [41]. We showed that the partition function of a 
general vectorlike theory can be split into a light and mirror part in an arbitrary gauge 
background, each part containing a chiral fermion representation. We also demonstrated, 
generalizing results of Neuberger [43] and Liischer [44, 45, 46] to arbitrary chiral actions, 
that while the partition function of the vectorlike theory is a smooth function of the gauge 
background, the splitting of the partition function is singular — its light and mirror parts each 
have singularities in gauge-field space — iff the mirror and light fermion representations are 
not separately anomaly- free. While these results indicated that anomalous and anomaly- 
free mirror theories differ in an essential way, as Golterman and Shamir pointed out soon 
thereafter [47], the singularity of the light-mirror split of the partition function does not 
resolve the puzzle alluded to above (see also the Addendum of [41]). This is because the 
singularity of the light-mirror split is topological and thus there is freedom — not unlike the 
ability to move the location of a Dirac string — to render the split of the partition function 

*We will liberally use this term in the two-dimensional context here, instead of the more appropriate but 
unwieldy "lack of algebraic order." 

^For example, a nonunitary theory, where the imaginary part of the Euclidean polarization operator is 
nonlocal, but the real part is local, see Section 3.3. 
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locally smooth at any point in gauge field space. Hence, as will be further expanded upon in 
this paper, there is still a puzzle in the case where the gauge field is non-dynamical, as in the 
simulations of [38]. 

As we will see in this paper, our previous results [41] provide us with important tools 
to further study the attempted decoupling of an anomalous representation. Since the two- 
dimensional model studied in [38] is the simplest theory with nontrivial mirror dynamics and 
is relatively inexpensive to simulate, understanding the answer to our question in detail is 
bound to give us clues as to the working of the EP proposal formulated with exactly chiral 
lattice fermions. In particular, we hope to distinguish between the possibilities (a. and b.) 
outlined two paragraphs above. 

1.4 Outline and summary 

We begin in Section 2, by first describing the deformation of the Schwinger model, studied in 
[38], that we call the "1-0" model. While this model is the subject of our numerical studies 
here, the analytical results of this paper are more generally applicable. In the end of Section 2, 
we review the main results of [41], notably the "splitting theorem" about arbitrary variations 
of general chiral partition functions with the gauge background, as this technical result is one 
of the major tools making this study possible. 

The considerations of the following Section 3 are quite general and apply to any EP- 
like model in either two or four dimensions. We first show that gauge invariance of the 
full partition function, combined with the local smoothness of the light-mirror split, implies 
certain anomaly matching conditions on correlation functions of gauge currents in the mirror 
theory. In two dimensions, the simplest nontrivial condition is on the two-point function 
in vanishing gauge background, i.e. on the polarization operator.^ The anomaly-matching 
conditions on two- or three-point functions of the mirror gauge current in vanishing gauge 
background are exact and do not depend on the details (i.e. strength) of the mirror-sector 
interactions — a simple consequence of the "splitting theorem" [41] . We note that the anomaly 
matching conditions do not require dynamical gauge fields and introducing infinitesimal gauge 
backgrounds is sufficient to establish them. 

The conditions derived in Section 3.3 are the equivalent of 't Hooft anomaly matching 
in theories with strong infrared dynamics: whatever the strong mirror dynamics, the mirror 
spectrum — assuming long-distance unitarity — has to be such that these anomaly matching 
conditions hold. A further similarity is that these conditions can be also usefully applied to 
symmetries that are not gauged in the target theory.'' The fact that conditions like 't Hooft 
anomaly matching should be obeyed by a strongly-coupled mirror theory on the lattice should 
not come as a surprise — that they can be precisely formulated and studied already in finite 

^In four dimensions with zero gauge background one would have to consider the three-point function to 
exhibit the effect of the anomaly; however a study of the polarization operator would still be of use as a probe 
for the existence of massless charged states. 

^In Section 1.5, we give an example showing how anomaly matching implies that a seemingly viable mirror 
theory for the anomaly-free "345" chiral U{1) gauge theory is bound to produce massless mirror modes. 
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volume and lattice spacing is an immediate consequence of the existence of exact chiral lattice 
symmetries. 

As usual in strongly-coupled theories, finding out how anomaly matching conditions are 
obeyed is a nontrivial dynamical question, to which we turn next. In Section 3.3, we review the 
Goldstone and massless chiral fermion modes of solving anomaly matching and explain how 
we expect to distinguish between the two in our numerical simulations of the 1-0 model. Then, 
in Section 3.4, we list a number of exact properties of the polarization operators of arbitrary 
chiral theories (such as our mirror theory), which are independent of the couplings and are 
derived in Appendix A. The verification of these properties in a Monte Carlo simulation is 
used to provide important checks on its consistency. 

In order to learn how anomaly matching is realized in the 1-0 model, in Section 4 we begin 
the study of the mirror polarization operator. We first present the rather technical derivation 
of the mirror-theory gauge two-point function. Since the details hold for general theories and 
may be of interest for future studies, we explain them in several steps in Sections 4.1 and 
4.1.1. The final expression of the mirror polarization operator in terms of correlation functions 
of the mirror theory (to be computed via Monte Carlo methods) is given in Section 4.1.2. 
Further details needed to express the correlator via variables used in the actual simulation 
are given in Appendices B, C, D. 

As already stated, our main interest is in the polarization operator of the mirror theory. 
The real part of the polarization operator at long distances contains information on the exis- 
tence and number of massless charged degrees of freedom of the mirror theory. In Euclidean 
space, the imaginary part contains the mirror contribution to the anomaly. As a check on our 
simulation, we have, in each case listed below, verified that the divergence of the imaginary 
part is exactly as required by anomaly matching (and is thus equal to the negative of the 
divergence of the free light fermion polarization operator). 

The Monte Carlo results of this paper are presented in Section 4.2. We begin by first 
explaining our strategy for looking for massless poles in the polarization operator on an 8 x 8 
lattice^ via its small-momentum discontinuity as a function of direction. In Fig. 1, we show 
the small-momentum discontinuity of the real part of the light-theory polarization operator — 
that of a free chiral GW fermion, defined in eqn. (A. 5) of Appendix A. In later Sections, we 
compare the analytic free-fermion result of Fig. 1 to the Monte Carlo results for the strongly 
interacting mirror sector. 

In Section 4.2.2, we show the results of the Monte Carlo study of the mirror theory in 
the strong-coupling symmetric phase. This phase — where all mirror-sector global symmetries, 
apart from the (would-be) gauge symmetry, are explicitly broken by the non-gauge mirror 
interactions — would be of most interest for mirror-decoupling a la EP in anomaly-free models. 
We compute the mirror polarization operator and look for the presence of small-momentum 
discontinuities. The results for the real part of the mirror polarization operator for different 
values of the mirror couplings, shown in Figs. 2, 3, 4, and 5, are strikingly similar to the 

*As we explain later, the computational demands of the problem limit us to a rather small lattice. 
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free-fermion result of Fig. 1. The numerical values and small- momentum discontinuity of 
the mirror polarization operator indicate that the number of massless charged degrees of 
freedom in the strong-coupling symmetric phase is the minimal one needed to satisfy anomaly 
matching. These results show that the strong mirror dynamics obeys 't Hooft anomaly 
matching and, in this phase, gives rise to a single massless charged chiral fermion, realizing 
the possibility b.) of Section 1.3 via the minimal solution with massless fermions. 

In Section 4.2.3, we study the mirror theory in the broken phase, where the "spin-spin" 
coupling K is large, see (2.2). The discontinuity of the real part of the polarization operator 
in Fig. 6 is now consistently interpreted as due to a massless scalar ("Goldstone boson") 
and shows that 't Hooft anomaly matching is obeyed in the Goldstone mode. This is easily 
explained, as the physics in this phase can be understood in perturbation theory, even at 
strong Yukawa coupling [38] . 

In Section 4.2.4, we study the strong-coupling symmetric phase with vanishing Majorana 
Yukawa coupling, /i — > of eqn. (2.1). The real part of the polarization operator, shown on 
Fig. 7, has, as in the h > 1 phase, a discontinuity like that of the free chiral fermion. However, 
the numerical values of the real part of the polarization operator are approximately (we are 
in a rather small volume) three times larger. This difference is due to the presence of three 
massless charged chiral degrees of freedom — a massless chiral fermion and a massless vectorlike 
pair.^ Thus, the spectrum in this phase is consistent with 't Hooft anomaly matching, but 
has more massless charged fermions than required by anomaly matching. This is explained in 
Section 5, where we also show that only at /i = we have analytic control over the symmetric 
phase. 

The final numerical result of this paper is presented in Section 4.2.5, where we study the 
scaling of the real part of the K-dependent contributions to the polarization operator deep 
in the strong-coupling symmetric phase (for decreasing values of k — > 0). As Fig. 8 shows, 
these terms do not contribute to the discontinuity of the polarization operator. Thus, the 
discontinuities, and hence the massless fermions, in the strong-coupling symmetric phases 
arise from the terms containing solely mirror-fermion current-current correlators. 

In Section 5, we use another representation of the partition function, obtained via a field 
redefinition. This representation is particularly useful to analytically interpret the results 
from Section 4.2.4, in the vanishing Majorana coupling strong-coupling symmetric phase 
(unfortunately, this is not the phase which would be useful for the study of decoupling in 
anomaly- free models), and explains the origin of the three massless modes found numerically 
in Fig. 7. 

For convenience of the reader, our conclusions are given in Section 1.5 below. We stress 
again the main results of this paper — that 't Hooft anomaly matching holds in lattice theories 
with strong non-gauge dynamics and (assuming unitariry) imposes important constraints on 
the spectrum at strong coupling. Our Monte Carlo study of the "1-0" model mirror dynamics 
shows that the anomaly matching conditions are obeyed by the mirror theory. In every 

®The massless vectorlike pair does not contribute to the anomaly; as a consistency check, we have verified 
that the divergence of the imaginary part of the polarization operator is as required by anomaly matching. 
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region of the phase diagram we studied, the strong GW-fermion mirror dynamics gives rise 
to massless charged states, reahzing option b.) of Section 1.3 — in a manner consistent with 
unitarity of the long-distance theory. Finally, we discuss the implications of these results for 
anomaly-free models and outline directions for future study. 

1.5 Conclusions, and where do we go next? 

The discussion of this paper shows that when a fermion formulation with exact lattice chirality 
is used, the question of decoupling the mirror fermions from a vectorlike gauge theory is 
intimately intertwined with their contributions to anomalies. While this conclusion is not 
unexpected, these issues have usually not been central to past studies (including our own) of 
attempts to decouple the mirror fermions in strong-coupling symmetric phases. Such studies 
were usually confined to purely non-gauge strong dynamics at the scale of the lattice spacing, 
without much emphasis on possible mirror-sector contributions to anomalies. This is largely 
because correlators of global chiral currents in the mirror theory could not be studied, as 
the exact lattice chiral symmetries (which enabled us to establish many exact properties of 
these correlators already at finite lattice spacing and volume) were not known. The main 
contribution of this paper is the precise formulation on the lattice of 't Hooft's anomaly 
matching conditions, applied to strong mirror interactions, and the study of their solution in 
a particular example. 

This paper shows that models with mirror Yukawa interactions formulated via exact 
lattice chirality are "smart" enough that, consistent with anomaly matching, in the strong- 
coupling symmetric phase they lead to massless degrees of freedom, rather than to a nonuni- 
tary theory. In our simple toy model, the minimal number of chiral fermions needed to cancel 
the light fermion anomaly remain massless in the strong-coupling symmetric phase. The mir- 
ror spectrum emerging from the cutoff-scale strong dynamics obeys 't Hooft anomaly matching 
and is consistent with the preservation of (would-be-gauged) unbroken global symmetries. 

A final obvious question is what our results here imply for anomaly-free models. Consider, 
for example, the "3-4-5" model, a two dimensional chiral U{1) gauge theory with 3_, 4_ and 
5+ massless fermions (the number indicates the U{1) charge and it the chirality). A direct 
argument predicting massless mirror states, based on the nonvanishing gauge anomaly of the 
light sector, can not be applied now. What our current knowledge allows us to say for sure 
is that whether there are massless mirror states — or not — depends on the implementation of 
the strong mirror interactions, in particular on their symmetries. 

The anomaly matching conditions find another use in this regard, as we now explain. 
Consider implementing the decoupling or mirrors from a vectorlike theory with three Dirac 
fermions of charges 3, 4, 5 by adding three uncharged mirror fermions and three unitary 
scalars, i.e. taking three copies of the 1-0 model and appropriately changing charges/chiralities 
(one could call the three copies the 3-0, 4-0, and 5-0 models). From our comments above, it 
is clear that such an implementation of the 345 model would lead to massless mirror states 

^°Thus, comparing Fig. 1 to Figs. 2-7 can be considered a Monte Carlo "proof of 't Hooft anomaly matching. 
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already in trivial gauge background, despite two facts that might suggest otherwise: that the 
mirror and light spectra are anomaly- free and that, when the global U{1) appropriate to yield 
the 345 model is gauged, both the anomaly-free and anomalous global symmetries of such a 
lattice implementation are as in the target continuum theory. The easiest way to see that 
massless modes will result at zero gauge background is to weakly gauge the three separate 
global U{1) symmetries present in the 3-0, 4-0, and 5-0 models and demand consistency of 
the two-point current correlators, as we did in this paper. 

A general way to phrase the condition that the mirror interactions should obey in order to 
avoid massless states due to anomaly matching from an extra global symmetry is to demand 
that when the gauge interactions are turned off, the mirror theory should have no global 
symmetries other than the global part of the gauge group. With gauge interactions turned 
off, the implementation of the 345 model of the previous paragraph has two extra global 
C/(l)s, which act simultaneously on the light and mirror components and thus impose further 
conditions on the mirror spectrum implying the existence of massless modes. 

This example leads us to conjecture that if the mirror interactions couple the 3+, 4+, 
and 5_ mirrors by adding only one scalar and a single neutral 0_ mirror fermion (needed to 
have a sensible static limit, see [23]), including the most general gauge- invariant couplings 
breaking all mirror global chiral symmetries, there wouldn't be any massless mirror states 
in the strong-coupling symmetric phase. At the moment, our strongest argument is that 
with all the mirror global symmetries explicitly broken, there is no reason we know of, such 
as anomaly cancellation of any symmetry, for massless mirrors to exist for all values of the 
couplings. Conversely, if a future study finds that massless mirror fermions exist for all values 
of the mirror couplings, absent any global chiral symmetry in the mirror sector, we expect 
that there should be some other reason for this — a symmetry or a deeper dynamical principle 
we are not yet aware of. Needless to say, it would be of great interest to understand what 
this principle might be. 

Finding out whether our conjecture above is true is left for future work. Two avenues for 
progress seem promising: i.) coming up with a theoretical argument why decoupling should 
always fail and ii.) verifying or refuting our conjecture via a numerical "experiment", similar 
to the one of this paper, but this time with an anomaly-free model. Since no dynamical gauge 
fields are needed at this stage, the study of (say) the 345-mirror dynamics is quite feasible 
given appropriate computer resources and the techniques already developed in [38], [41], and 
the present paper. 

2. The "1-0" model and the "sphtting theorem" 

The Yukawa-Higgs-GW-fermion model considered here, which we call the "1-0" model, is 
a U{1) two-dimensional lattice gauge theory with one Dirac fermion of charge 1 and a 
neutral spectator Dirac fermion x- Considering this theory is motivated by its simplicity: it 
is the minimal Higgs- Yukawa- GW- fermion model in two dimensions which holds the promise 
to yield, at strong Yukawa coupling, a chiral spectrum of charged fermions and is, at the 
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same time, amenable to numerical simulations not requiring the use of extensive computing 
resources. The fermion part of the action of the "1-0" model is: 

S — Slight ~l~ Sqjiirror (^-l) 

Slight = - • -Di • -0+) - {X- ■ Do ■ X-) 

Smirror = " (V'- ■ Di-%l)J) - {x+ ' Dq ■ X+) 

+ y • 0* • x+) + (x+ • • V'-) + /i ■ ^12 ■ x+) - {x+ ■i2-<f ■ } . 

The chirality components for the charged and neutral fermions are defined by projectors 
including the appropriate Neuberger-Dirac operators^^ (charged Di and neutral Dq) for the 
barred components, i.e. •i/'± = V'(l -F Ts)/^ • The brackets denote summation over the lattice 
sites as well as a spinor inner product. The field (j)x = e^^'^ , \'qx\ < vr, is a unitary higgs field 
of unit charge with the usual kinetic term: 

= f E E [2 - ( '^x Ux,.+(. <A.+A + h.c. )] . (2.2) 

The inclusion of both Majorana and Dirac gauge invariant Yukawa terms is due to the re- 
quirement that all global symmetries (including those of the mirror fermions) not present 
in the desired target chiral gauge theory be explicitly broken, see [23]. Moreover, consistent 
with the symmetries, if the Majorana coupling h vanishes, there are exact mirror-fermion zero 
modes for arbitrary backgrounds , which can not be lifted in the disordered phase [38] . The 
lattice action (2.1) completely defines the theory via a path integral over the charged and 
neutral fermion fields, the unitary higgs field, as well as the gauge fields. We will not perform 
the integral over the lattice gauge fields, but will study in detail the variation of the partition 
function with respect to the gauge background. 

From now on, we will call the fermion fields that participate in the Yukawa interactions 
the "mirror" fields — these are the negative chirality component, of the charged ip, and 
the positive chirality component, x+i of the neutral x — while the fields ■;/'+ and X- will be 
termed "light." 

Our interest is in the symmetric phase of the unitary higgs theory, expected to occur at 
K < Kc — 1, where the higgs field acts essentially as a random variable (modulo correlations 
induced by k 7^ and by fermion backreaction) . Based on experience with strong- Yukawa 
expansions in theories with naive or Wilson fermions, it is expected that in the large-y, &xed-h 
limit, there is a symmetric phase. 

The analysis of the 1-0 model of ref. [38] was performed in a vanishing gauge background. 
The eigenvectors of 75 were explicitly worked out (see also Appendix C) and used to manifestly 
split the partition function into "light" and "mirror." A Monte Carlo simulation of the mirror 
sector at infinite-y and &xed-h was performed. For small values k < Kc a strong-coupling 



^See Appendices B and C for our convention. 
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symmetric phase exists for h > 1 and for h ^ 0, while at h ~ 0.7 evidence for a Berezinksii- 
Kosterhtz-Thouless-hke transition was found. The spectrum of the mirror theory was also 
numerically probed by studying correlators of local operators. The evidence found — using 
local operators to probe the charged fermion spectrum — pointed toward decoupling of the 
mirror sector, with no breaking of the chiral symmetry of the mirror sector (the symmetry to 
be gauged by the U{1) gauge field). 

These results led to the puzzle already alluded to and to the studies of [41] and the 
present paper. In [41], we studied the splitting of the partition function into "light" and 
"mirror" parts also in an arbitrary nonvanishing gauge background. To remind the reader 
of our convention (identical to that used in the simulation of [38] but slightly different from 
the one of [41]), we use the definite-chirality eigenvectors of 75 and the projectors P± on the 
corresponding spaces: 

%Ui = -Ui , %Wi = Wi , (2.3) 
i i 

where we treat u,w as columns and u^,w^ as rows. We also use the eigenvectors of 75 (the 
latter are independent of the gauge background) and the associated projectors P±: 

l5Vi = Vi , 75ti = -ti , (2.5) 
P+ = Y,v^vl , P. = Y,t^t\ = l-P+ . (2.6) 

i i 

Using (2.3), (2.5), a general Dirac field ^'x, can be decomposed into chiral components af , 
via the 75 eigenvectors, while the conjugate spinor field "ifx is decomposed into chiral compo- 
nents af,a^ using the 75 eigenvectors, as follows: 

^x = J2 + , (2.7) 

i 

= ^aV'ul(x) +aLu;l(x) . (2.8) 

i 

Applying the split (2.7) to the fields ifj and x of the "1-0" model, we note that only the 
charged eigenvectors (of both light and mirror fields) depend on the gauge background. The 
expansions (2.7) of the "mirror" fields are explicitly given below: 

X+ = ^/J>i, x+ = E^W[0], (2.9) 

i i 
i i 

^■^ Analytic [48] and Monte Carlo [49] studies of a similar four-dimensional theory with no Majorana coupling, 
i.e. at h — 0, also found evidence for a strong-y symmetric phase. 
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Clearly, expansions similar to (2.9) hold for the "light" fields as well: 

X- = Y^PlU , X- = Y.^-wl[U] , (2.10) 

i i 
i i 

After substitution of (2.9, 2.10), the partition function of the model (2.1) splits as follows: 
Z[U-yM = Z+[U]^ Z^[U-yM . (2.11) 

Here Zj^\U\ = det ||(?iJ[J7] • D\U\ ■ f j ) 1 1 x (similar determinant for the neutral light spectator 
X-) is the light sector partition function. The jacobian J is a product of jacobians for the 
charged and neutral sectors; see [41] for details. The mirror partition function is denoted by 
Z_ and is given, more explicitly, by an integral over the charged mirrors (a_,a_), neutral 
mirrors (/3+,/3, ), and unitary scalar field: 



Z.[U;y,h] = j (fa- (f(3+ e"^-""'"°'-^« , (2.12) 



where the mirror action from (2.1) is expressed in terms of the integration variables a_,/3+ 
and the eigenvectors via (2.9). The mirror fermion integral is thus a determinant which 
includes the kinetic term and Yukawa terms from (2.1) and the mirror partition function is 
the average of the determinant over the random (in the disordered k — > phase) unitary field 
(f>x (in (2.12), d<j) denotes a path integral over the phases of (p). The mirror partition function 
(2.12) with U = 1 was the object used to calculate local mirror observables in the simulations 
of ref. [38]. 

When U ^ 1, the mirror partition function Z_ (2.12) depends on the gauge background 
through the operators entering Smirror + S^, (the Neuberger-Dirac operator and the associated 
projectors that appear in (2.1)) as well as through the gauge background dependence of the 
eigenvectors of 75 used to split the partition function (zi;j[^], see (2.9)). This dependence 
was studied in [41], where an important technical result was derived: under an arbitrary 
variation of the gauge background, the variation of the mirror partition function, due to the 
variations of both the eigenvectors and the operators entering the action, factorizes no matter 
how complicated the mirror partition function is. Explicitly, the "splitting theorem" states 
that for an arbitrary variation of the gauge background: 

6logZ.[U] = J2iSwlw,) + 6o\ , (2.13) 

i ^ ' 

where "( • )" denote an expectation value calculated with the partition function Z_ and O 
collectively denotes the various operators depending on the background which appear in 
the action (the Neuberger-Dirac operator D and corresponding projectors P±). This is an 
important result, as it encodes on the lattice the idea that anomalies do not depend on the 
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action (see the discussion in [41]). Furthermore, as we will see later, the splitting theorem is 
indispensible in the calculation of mirror gauge-current correlators in a perturbative expansion 
in the gauge field. 

Coming back to the full 1-0 model partition function, we note that because the l.h.s. of 
(2.11) is manifestly gauge invariant, so is the r.h.s., since it is obtained from the l.h.s. simply 
via a (locally) nonsingular change of variables. We know how two of the factors on the r.h.s. 
transform under gauge transformations — the light partition function Z^[U] and the Jacobian 
J[C/]~^, for which we have from [41], for infinitesimal gauge transformations uj (as indicated 
by the ~ sign): 



where d^jUi denotes the variation of i-th eigenvector of 75 under a gauge transformation. 
Therefore, from (2.14), and the fact that the l.h.s. of (2.11) is gauge invariant, it follows that 
the mirror partition function transforms under gauge transformations as: 



independent not only on the values of the Yukawa couplings (y, h) but also most of the details 
of the mirror action. We also note that the "splitting theorem" (2.13), valid for general 
variations, gives a direct proof of (2.15) when restricted to gauge transformations [41]. 

The gauge variation of the mirror partition function leads us to the already mentioned 
paradox. The exact result (2.15) shows that the gauge transformation of the mirror partition 
function is independent of the Yukawa coupling and should precisely cancel that of the light 
chiral fermion (2.14). If, at 1/ ^ 00 and h > 1, the mirror sector only involves heavy degrees 
of freedom, as the numerical results of [38] suggest, and if these zero-background results 
persist for arbitrarily small gauge backgrounds (as one is inclined to expect based on local 
smoothness of the basis eigenvectors), then the mirror partition function should be a local 
functional of the gauge background. But eqn. (2.15) argues that this local functional's gauge 
variation must precisely cancel the anomaly of the light chiral fermion, which is known to be 
impossible as the anomaly is not the variation of a local functional. 

This question was a major motivation for the study of [41], where we showed that the 
"light" -"mirror" split of the partition function is a singular function of the gauge background 
in any model where the mirror matter representation is anomalous. In particular, the mirror 
partition function (2.12), used in the numerical simulation in [38], has a discontinuity precisely 
atU = 1. We also conjectured there that the singularity of the mirror partition function might 
play a role in resolving the paradox, as the results on the mirror spectrum obtained in [38] 
via the singular Z_ (2.12) may not survive an integration over the gauge field or even a 
perturbative expansion around U = 1. While this claim may be plausible, if the gauge field 
is considered as an external background, with small fluctuations around it taken into account 




(2.14) 




(2.15) 
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perturbatively, a paradox still persists, as is briefly explained in the Addendum of [41]. This 
is the issue we want to address here. 



3. Smoothness, the light-mirror split, and anomaly matching 

We begin this Section by stating in a more formal way the conflict between the numerical 
results of [38], which found no evidence of long-range correlations in the mirror sector at 
y — > cx) and h > 1 (when probed with local charged operators), and eqn. (2.15), which states 
that the mirror partition function of the 1-0 model should effectively act as a Green-Schwarz 
term canceling the light sector anomaly. Our goal here is to formulate the anomaly matching 
conditions in a way useful for further study. 

3.1 Transversality of the full partition function 

The partition function (2.11) of the 1-0- model is gauge invariant, i.e: 

In Z[A + 6^A] = In Z[A] , (3.1) 

where U{x,x + ^) = e*^''^^^ with 5^^^(x) = — V^w^;, and V ^oJx = ^x+^l — ^x- This implies 
that: 

fix ' 

where y],iOx = lOx - i^x~fi- Taking of (3.2) gives: 

VV* . ^ . / , =0 , 3.3 

which, by expanding in Ap, around Ap^ = 0, implies transversality of all n-point functions: 

b'^ZXA] 



bAp{x)bAp^ (xi) . . . bAp^_^ {xn-\) 



= . (3.4) 

^=0 



Eqn. (3.4) should apply to the full partition function of the 1-0 model (2.1). The derivation 
of it assumed that Z\A\ is a smooth function of the gauge potential in the vicinity of vl = 0. 
Smoothness of Z\Ai\ holds because the vectorlike theory has a well defined measure and 
an action which is smooth with respect to the gauge background. (We note again that the 
singularity discussed in [41] appears separately in the "light" and "mirror" partition functions; 
the total partition function is that of the anomaly- free vectorlike theory and is nonsingular.) 



^^We stress again that smoothness holds subject to the "admissibihty condition" on the gauge background 
[36, 37]. 
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3.2 Local smoothness of the hght-mirror split 

The singularity in the "light" -"mirror" split of the partition funciton [41] is of topological 
nature. In the case when the "light" and "mirror" sectors are anomalous, there is a topological 
obstruction to defining their separate fermion measures (or of the "measure current" of the 
75 eigenvectors) as a globally smooth function of the gauge background. The mirror partition 
function and correlators studied in [38] depend on the 75 eigenvectors, which are discontinuous 
functions of the gauge background when it is turned on. 

The topological nature of the singularity means that its location in the space of gauge 
backgrounds can be moved around by redefining the phases of the basis vectors (equivalently, 
of the measure current). An explicit example is discussed in [41] for the homogeneous Wilson- 
line subspace. More generally, due to the smoothness of the P projectors, a locally smooth 
basis of eigenvectors exists at any given point of the space of gauge backgrounds. The possi- 
bility of choosing a basis which is locally smooth implies that in an expansion around Afj_ = 0, 
(3.4) can be applied to the split partition function, as we discuss explicitly below. 

In what follows, we present a lengthy derivation of the anomaly matching condition. We 
will use the results obtained here to calculate the contributions to the polarization operator 
of the gauge field from the "light" and "mirror" fields in later sections. To begin, we define 
the polarization operator of the full theory as: 



11^1^(2;, y) 



6HnZ[A] 



(3.5) 

A=0 



5Af,{x)6A4y) 

Since the full partition function splits as in eqn. (2.11), we have that: 

In Z[A] = In Z+[A] - In J[A] + In Z_ [A] . (3.6) 

Since the split is locally smooth in the neighborhood of A^ = 0, eqn. (3.6) defines a split of 
the polarization operator (3.5) into "light" (H"*"/"^, with the Jacobian contribution included 
into the "light" sector) and "mirror" (n~) parts: 

n^,(x, y) = U^l^ix, y) + H-,{x, y) . (3.7) 

Consider first the "light" polarization operator. Under an arbitary infinitesimal change 
6rj of the gauge background, the "light" partition function Zj^[A] = det{ul[A] ■ D[A] ■ vj) 
transforms as: 

In Z+[A] = Tr{P+D-^Sr,D) + J^(<5^nJ • uj) 

j 

= TT{P+D~'6rjD)+j;[A] , (3.8) 

where j^[A] is implicitly defined above. Similar formulae are derived in [41] (we only give 
this one in detail due to the different convention of this paper — here we use Pzf to define ip±). 
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We use "Tr" to denote trace over both spinor and space-time indices, and "tr" to denote 
summing over spinor indices only. We have also introduced the "measure currents" j^'"': 



(3.9) 



in terms of which the variation of the Jacobian is given as in [41] (modulo the change of 
convention, see eqn. (2.16) there): 



5,lnJ[^]=j-[A]+j-[A] 



(3.10) 



By combining (3.8) and (3.10), we find that the "light" plus Jacobian contribution to the 
change of Z[A] of (3.6) under a gauge transformation (cf. eqn. (2.17) of [41]) is: 



5^ In ^ = -j:[A] .M\ , 



(3.11) 



where j^j denotes the measure current (3.9) now corresponding to a gauge variation of the 
background. Now we take of (3.11) to find: 



ouJx 2 



Using the identity: 



J[A] 



eqn. (3.11) is clearly the same as: 

6\n[Z+[A]J-^[A]) 



8A^{x) 



5A^{x) 



(3.12) 



(3.13) 



(3.14) 



where we introduced the shorthand notation: 



6A^{x) 



(3.15) 



for derivatives 6^ to be used further (note that for brevity we often suppress the space-time 
index which we understand to be included in fi). Finally, we expand (3.14) around A^ = 
to linear order, allowed by local smoothness: 



6^lii{Z+[A]J~^[A]) 



5A^{x)5A^{y) 



-V 



A=0 



SA^y) 



_ i 6trj5xx[A] 
2 6A,{y) 



(3.16) 



A=0 
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Eqn. (3.16), using our definition (3.5, 3.7) of the polarization operator, is equivalent to: 



.14 



1 

A=0 ^ 



(3.17) 

^=0 



showing that the "light" polarization operator is not transverse. 

There are two contributions to the r.h.s. of (3.17): the first term, proportional to the 
derivative of the "measure current," exactly cancels with the identical contributions of the 
"mirror" sector — see eqn. (4.10) for the mirror polarization operator. The second term on the 
r.h.s. of (3.17), proportional to the derivative of the topological lattice field tr752;2;, represents 
the anomaly of the "light" fermions.^^ To make contact with the anomaly in the continuum, 
we note that the topological lattice field can be expressed as (the four-dimensional proof of 
[42] is trivially downgraded to two dimensions): 

tr 75.. = -^e^uF^'" + V*/i^[^] , (3.18) 

where F^j^i, = V ^Ay{x) — V yA^{x) is the field strength of the gauge potential A^{x) = 
—ilnU{x, fi), h^^[A] is a gauge invariant local current, and ei2 = 1. An explicit form of 
/i'^[74] can be obtained with some work; for example, the part of h linear in A and valid for all 
momenta, can be derived using the formulae in Appendix D, which also presents a derivation 
of the first term in (3.18). 

3.3 Anomaly matching and its possible solutions 

By the local smoothness of the "light"- "mirror" split and by gauge invariance of the full parti- 
tion function, the divergence of the "mirror" polarization operator should exactly cancel 

(3.17) . Since the measure current dependent parts of the "light" and "mirror" polarization 
operators always cancel, the important contribution to the divergence is the one containing 

(3.18) and representing the true effect of the anomaly. We write the cancellation requirement 
in the form — taking the low-momentum limit and implying a sum over repeated indices: 

iqf'fl-M = ^^uxq^ + 0{q^) , (3.19) 

where IT denotes the appropriately defined Fourier transform of the polarization operator — 
see (D.l, D.3, D.5). Now of course the rhs of (3.18, 3.19) is local, and one wonders if the 
usual continuum argument that it can not be the divergence of a local expression applies on 

^''We include the Jacobian contribution into the light polarization operator; note that in accordance with 
the "splitting theorem" it will cancel with a similar contribution from the polarization operator of the mirror 
partition function, see [41]. 

While equations similar to (3.17) will hold for the non-transverse higher derivatives of the light partition 
function, considering only the polarization operator in a trivial gauge background is sufficient to study the 
interplay between the anomaly and the light degrees of freedom in 2d. In 4d the trivial gauge background 
analysis would have to be extended to the three point function in order to capture the effect of the anomaly. 
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the lattice. A quick argument showing that it does is as follows. Rewrite (3.19) as the set 
of two equations for the imaginary part of the polarization operator, denoted with the same 
symbol for brevity, with c = — i/(27r): 

n- = _(c + fl-) ^ , (3.20) 

nri = (c-n2i)^ = (c-nr2) 

where we used II12 = due to local smoothness. Locality of XIn and 1122 would then 
require that c — II12 = ^Qi + c + 11^2 = ^Q2 + where A and B are arbitrary constants 
and dots denote higher powers of momenta, leading to: 

n- = -c + 5^2 + 0(9^) , (3.21) 
fli2 = c- Aqi + 0((?2) , 

conditions which are clearly incompatible. 

On the other hand, a nonlocal solution of the anomaly conditions (3.20) for the imaginary 
part of the polarization operator is given by = ^{Q2 ~ 9i)/(9i + li)- This is, in fact, the 
continuum value of the non-transverse part of the polarization operator for an anomalous 
theory, see eqns. (3.23, 3.24) below. A local solution of (3.20) can be found only if ni2 / 1121, 
i.e. only if there is a singularity near = such that the second derivatives of Z_ do 
not commute; this, however, goes against local smoothness of the measure current. Thus, 
our conclusion is that local smoothness combined with the nonvanishing anomaly of the light 
fermion (3.17) imply that the imaginary part of the "mirror" should have a nonlocal 
part. Eqn. (3.19), which leads to this conclusion is a mathematical consequence of the gauge 
invariance of the 1-0 model partition function and the freedom to choose locally smooth 75 
eigenvectors. 

We now enumerate the possible solutions of (3.19). First, note that in Euclidean space 
the anomaly appears in the imaginary part of the polarization operator and this remains true 
on the lattice, for arbitrary chiral theories formulated with GW fermions (see Appendices 
A, D). Next, it is well-known from the continuum that in a unitary Lorentz (Euclidean) 
invariant theory the zero-momentum singularity in the solution of (3.19) is due to either a 
massless Goldstone boson or a massless fermion [50, 51]. The massless scalar or fermion will, 
of course, also give a nonlocal contribution to the real part of the polarization operator, in 
addition to the nonlocal contribution to the imaginary part required by the anomaly. On the 
other hand, if our complex Euclidean partition function had no Hamiltonian interpretation, 
i.e. led to a non-unitary long distance theory, one could imagine that the imaginary part has 
a nonlocal contribution while the real part does not. 

Our strategy to look for massless charged mirror modes will be to study both the real 
and imaginary parts of the polarization operator of the mirror theory. We already know 
that its imaginary part has a nonlocal contribution giving rise to (3.19). The real part of 
the polarization operator probes the number and nature of the massless charged degrees 
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of freedom and thus gives information of the spectrum. In particular, in two dimensions, 
massless (Green-Schwarz) scalar or fermion loops lead to = poles in the real part of the 
polarization operator, see (3.23, 3.24). It is this part of the mirror polarization operator that 
will be of most interest to us. 

To recall how this plays out in the continuum and see what we might expect from our sim- 
ulations, consider first the example of a Green-Schwarz scalar theory with Euclidean partition 
function: 

Zgs[A] = Jvr] e/^'^(-f(9A'^-^M)'+i^^i^i2) ^ (3.22) 

such that the partition function is not invariant under gauge transformations. The normal- 
ization of the kinetic term is chosen such that it is the naive continuum limit of of (2.2), 
with (f> = e*^; note that perturbation theory is good when k ^ 1. We define the gauge boson 
polarization operator as in (3.5), and explicitly compute its Fourier transform from (3.22), 
with the result: 

The polarization operator is not transverse and its divergence is iQn^'Gs\A=o^'^^ ~ ^^^''^p^ 
in accordance with (3.19). The nonlocal imaginary part of (3.23) is precisely equal to the 
solution found in the paragraph below eqn. (3.21). The real part of (3.23) is the contribution 
of the massless scalar to the gauge field effective action; notice the shift of the coefficient 
K — > K -|- -^^T^ due to the anomalous Green-Schwarz term. If our mirror theory has a massless 
Green-Schwarz scalar, it would have to manifest itself by contributing to both the real and 
imaginary parts of the mirror polarization operator, as in (3.23). In fact, we will find (not 
unexpectedly) that at large k ^ 1, i.e. in the "broken" phase of the unitary Higgs mirror 
theory, this solution of anomaly matching is realized. 

A massless mirror charged chiral fermion would, similarly, also contribute to both the 
real and imaginary parts of II^i,: 

1 ( q^q" i e-'PqpqP + ePPqpq" 



nfL„U..(«) = TA-t- '"J - ^ '-^"^^^ . (3.24) 

where the real part of Ilpuiq) of the chiral fermion is equal to one-half that of the Dirac 
fermion in the Schwinger model. Note that, as opposed to the Green-Schwarz scalar (3.23), 
the coefficients of the real and imaginary parts of the massless chiral fermion 11^,^ are the 
same — their ratio is 1 vs. ~ 2'kk for the scalar when k ^ 1. Again, running ahead, we 



^^The regularization ambiguities discussed in [39] do not arise here, since (3.24) is the small-momentum 
limit of the contribution of a free chiral GW fermion to the basis-vector independent part of the polarization 
operator, see Appendices A, B, D. It should be clear that we do not claim that (3.24) represents the polarization 
operator of a consistent unitary anomalous theory ((3.24) would correspond to the singular a — 1 case of [39]). 
Instead, eqn. (3.24) is exactly what a massless mirror at y = would contribute to the full polarization operator 
of the vectorlike theory. 
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will find strong evidence that in the disordered phase at k < Kc the mirror spectrum obeys 
anomaly matching via massless chiral fermions. 

In our further study, we will compare our Monte Carlo results for the mirror Rell^jy with 
the lattice analogues of eqns. (3.23) and (3.24), which we can easily compute at finite volume 
and lattice spacing for a free charged unitary scalar and a GW fermion, respectively. We 
will see that this comparison can already be made on an 8 x 8 lattice, providing compelling 
evidence of what mode the mirror theory chooses to obey the 't Hooft conditions. 

3.4 Exact properties of chiral polarization operators 

We now list some exact properties of polarization operators that hold for general chiral the- 
ories, in particular for our mirror theory. These properties are derived in Appendix A. All 
equations below refer to the polarization operator in x-space and not to their Fourier trans- 
forms; recall that we absorb the space-time indices into fi, v (see (3.15)). 

The definition of 11"^ = 5^(^i,logZ_, see (3.5, 3.6, 3.7), implies that n~j^ is symmetric 
due to local smoothness of Furthermore, can be decomposed into a part that is the 
measure current derivative and the rest, see (4.7) and (4.10): 

n;. = 5,j;r + n;;. (3.25) 

As shown in Appendix A, n~j^' is always a total derivative: 

n;; = <5.n; ' . (3.26) 

and, in addition, n~ ' is exactly gauge invariant. Therefore, proceeding as in the derivation 
of eqns. (3.12, 3.14), we find: 

v:n;; = , (3.27) 

while with respect to the first index, we have from (A. 16) in Appendix A: 

v;n;; = UmiI,. (3.28) 

Now, the total mirror 11"^ is symmetric, but byj"^ and 11"^,' are separately not, but obey: 

(n;; - n;;) = -<5,j- + = t^, , (3.29) 

where is the curvature of the measure current, explicitly given in Appendix A, which is 
a known local functional of the gauge field whose divergence V*.7^^j/ gives half the anomaly. 
These results imply that the symmetric and antisymmetric parts of 11"^,' each contribute half 
of the anomalous divergence (3.28) (or (3.19)). 

In our numerical simulation we calculate n~j^', as the measure-current part is exactly the 
opposite that of the light theory. Since all properties listed in this Section hold independently 
of the mirror action, in particular of the strength of the mirror couplings, the verification of 
(3.27), (3.28), and (3.29), in a numerical simulation at strong mirror couplings provides an 
important check on its consistency. 
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4. Is anomaly matching satisfied? How? 



This Section serves two main purposes. First, in Section 4.1 we describe the analytical work 
required to find an expression for the mirror polarization operator in terms of correlation 
functions of the mirror theory, to be computed via Monte Carlo simulations. Second, in 
Section 4.2, we present the results of our simulations. 

4.1 Setting up the calculation 

To appreciate the technical details of the calculation of the mirror polarization operator, we 
begin by noting that if one varies the full theory partition function first, without using a 
locally smooth light- mirror split, and then substitutes = and the corresponding = 
basis vectors to calculate the polarization operator, one finds that the gauge current involves 
terms that mix "light" and "mirror" states. 

This point is already evident in considering the gauge current in the y = "light" plus 
"mirror" theory, which is given simply by (ip ■ 5^D ■ ip). Substituting the expansions (2.7) of 
the un-barred spinors in terms of the 75 eigenvectors u,t {ip = a'_ti + a\vi) and of the barred 
spinors in terms of the 75 eigenvectors w,u {ip = a^_w\ + a!|_ii|), one finds that: 



i.e. the light-mirror cross terms do not vanish. The a±a=p cross terms, upon insertion in 
a path integral over the fermions ip^ip-, contribute "light" -"mirror" loops. This fact would 
make the calculation of current-current correlators rather difficult, especially in the presence 
of interactions (at y 7^ 0) treating different chiralities differently. In particular, if a similar 
variation was done also in the Yukawa interactions, there would also be "light"- "mirror" cross 
terms, making the calculation rather inconvenient, as the "light" contribution is calculated 
"by hand" and the mirror — via Monte Carlo simulations. 

Using the locally smooth (in the infinitesimal neighborhood of = 0) basis vectors, 
together with the "splitting theorem" of [41], helps achieve a separation of the "light" and 
"mirror" contributions to correlation functions which makes the calculation manageable. The 
results for the full partition function do not depend on the basis used, and it is a great 
convenience to have "light" and "mirror" contributions to current correlators separated. Our 
goal is to discuss the correlators in the y = 00 mirror theory that contribute to 11"^,. We use 
the results of [41] to simplify the calculation. However, we start with a simple example — the 
calculation of the mirror polarization operator at y=0 — to illustrate the essential technical 
points. 

4.1.1 Warm-up: mirror polarization operator at y = 

This warm-up gives a better idea how the "splitting theorem" can be used. In the end, 
of course, adding the "light" theory piece just reproduces the polarization operator of the 



(4.1) 
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massless vectorlike theory at ?/ = (in a somewhat "twisted" left-right separated way, to be 
useful later). 

The y = "mirror" partition function — ignoring the neutral mirror, since it will not 
contribute to the polarization operator when y = — is: 



Z_ = y" da-da-e-^ , (4.2) 
S = -alcJ_{w\- P+- D -tj) 

and its first variation wrt A^[x) (below, we use to denote 'W^da'^_da^_) is: 

.„Z_ = /A_.-V_(*„u,t.P,.«.*,,e- (4.3) 



The first line above equals Z_ Xlil'^At^I ' ~ ^-JJH by the splitting theorem (2.13). 

We note that the second line in (4.3) represents a different chiral partition function (as, 
by the criterion of [41], it is a partition function whose variation wrt w"^ is orthogonal to n'l") 
and the splitting theorem can be used while calculating the second variation of this term. 
Finally, the third line in (4.3) vanishes, because D ■ tj = D ■ P^tj = P+ ■ D ■ tj, combined 
with (5^75 • 75 = —75 • (5^75, implies that after P+ is pushed through 5^P+ it becomes P„ and 
annihilates . We then have from the above comments: 



+ / d^a- 



(4.4) 



{w\- P+■5^D■tj)e 



^eyit, we calculate the second variation. A very important technical point is that the 
splitting theorem can be applied iff the variation of the partition function due to the change of 
the basis vector is orthogonal to the opposite chirality basis vector. This is the reason we keep 
the factor of P+ in the second line in (4.4) even though it might appear as a tautology here; 
the price to pay is that we have to account for the projector's variation when we calculate 
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the second variation of Z_ , as we will see below. Thus we find: 



KZ- X i 




(4.5) 



We used the splitting theorem on the third line above, wrt the chiral partition function given 
by the second line in (4.4). Armed with (4.4) and (4.5) we can now compute the y = 
"mirror" theory polarization operator, defined by (3.5) with Z — > Z_: 



where brackets "( • )" denote expectation values calculated with the mirror partition function. 
It is immediately seen (by recalling [41]), using (a'a-^) = \ \{wl. ■ D ■ ti)~^\\ij = (4 • ■ wj), 
that the polarization operator of the mirror can be written as: 



Many of the projectors above can be dropped due to the Ginsparg- Wilson relation, but were 
left in for comparison with (4.6). 

Now, it is straightforward to compute directly from the expression for the y = 
charged- mirror partition function Z_ = det \ \{wj ■ D ■ tj)\\ and see complete agreement with 
(4.7); we leave this as an exercise for the reader. The calculation in this section was done in 
such detail in order to see the agreement with the much faster direct calculation alluded to 
above and to emphasize that the term with the variation of the projector in (4.7) is crucial 
for this agreement — this projector was inserted to allow use of the splitting theorem when 
performing the second variation of the mirror partition function. 

Of course, for the case of interest y ^ 0, we do not have an expression for Z_ as simple as 
det \ \{wl ■ D - tj)]], so the best we can hope for is to cast the mirror n~j^ into a form similar to 
(4.6), which will give the mirror polarization operator in terms of mirror theory correlation 
functions. These can be computed via a Monte Carlo simulation at zero gauge background, 
by first computing explicitly the various functions of position and momentum appearing in 
(4.6), such as {wl ■ duP+ ■ S^D ■ tj), and then using the existing code of [38]. 



TT— C -W 

+ (^{aiaia'La'L) - (aLai) (a^al)) {wj ■ 5^,D ■ tj){wl ■ 6^D ■ tg) 



(4.6) 




d.j]-: + Tii6^5,D ■ P_ • D-' ■ P+) + Tii6,P+ ■ 6^D ■ P_ • D'^ ■ P+) 
Tr{5^D • P_ • D^^ ■ P+ ■ 5uD ■ P_ • D^^ ■ P+) . 



(4.7) 
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4.1.2 Mirror polarization operator at y / 



The 1-0 model action is given in (2.1). Substituting (2.9) into the mirror action, we find the 
terms in the mirror action that depend on the gauge field: 



Sn^irror = " ct.al {w\[A] ■ D[A] ■ tj) + y ct_(3{ {w\[A] ■ P+[A] ■ <t>*Vj) 

- yh %a'_ {u]-i2 ■ cP* ■ Pl[A] ■ w* [A]) + S,[A] (4.8) 
+ — independent , 



where is the kinetic action for </>. The mirror partition function is defined in (2.12). For 
the first variation of In Z_ [A] we find (brackets ( • ) now denote expectation values with the 
full mirror partition function (2.12)): 



+ ( alai {wl ■ P+ ■ 6^D ■ t,)) 

+ |((0*.<5^C/-0)+h.c.) (4.9) 

- y { aipi {wj ■ P+ ■ 6^P+ ■ 4>*v,)) 

- yh {alpi {u]j2 ■ <A* • S^P+ ■ P+ ■ w*) , 



where the first two lines on the r.h.s. appear just as in (4.4). Now we have to perform a 
second variation, carefully use the splitting theorem several times, and finally compute = 
(5^(5^1nZ_[^]|^_Q, watching for the many cancellations. This calculation straightforwardly 
follows the already established rules and after some tedious algebra we obtain the mirror 
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polarization operator: 



n 



— () i"" 

+ f {{((>* ■SuSf.U ■(!))+ h.c.) 

,^2 



+ 



4 

K 

2 L 



6^U ■((>)+ h.c.]) 



c 



y {w\ ■ S,{P+5^P+) ■ (P*v,)) 



(4.10) 



yh (aL/3t (^72 • ^ ■ 5.{6^pI • P|) • w*)) 



+ (/i ^ I/)] 

+ y^ {(alpiiwl • (^^P+ • ^*Vj) + /iaL:5'+(^72 • 0* • S^P^ • w^I) ) x 



C 



yK 
2 L 



([((/>* • <5/,[/ • (/>) + h.c] aipiiwl ■ 6^P+ ■ (j)*Vj) + haipliuU2 ■ • S^P^ ■ w*) ) 



with (...)'^ indicating the connected part. The first three hnes are identical to (4.6). 

The form of the mirror polarization operator (4.10) simplifies somewhat in the y = oo 
limit of our simulation when many terms vanish. At infinite Yukawa, the mirror theory 
conserves the number of + fermions minus the number of — fermions, implying that all 
correlators with an unequal number of + and — fermions vanish as y — > oo. Therefore, we 
have a simpler expression for the mirror polarization operator (4.10), which we write as a 
sum of several terms: 



n 



A"^ \ y—>oo 



A _|_ n' TT' = TT^ -I- TT'^^ -I- TT'' 



(4.11) 



where the fermion current-fermion current contribution U^i, is: 



HI 



y {alPi {wj ■ 5,{P+5^P+) ■ rvj)) - yh {al(3'^ {u]j2 ■ 0* • 6,{6^Pi • P|) • w*)) 



atp^+iwl • SuP+ • ^*vi) + ha'LP\{u]-f2 ■ (/)* • S^P:^ -wl))) 



X 
C 



(4.12) 



the scalar current-scalar current contribution II^j^ is: 



([(0* • 5^U ■ <P) + h.c] [{(/>* ■ 6,U ■ (P) + h.c.]) 



c 



(4.13) 
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and, finally, the mixed fermion current-scalar current contribution IIuv is: 



= 

2 

+ ly)} . (4.14) 



{([(0* • S^U • 0) + h.c] [aipiiwl • • (l)*vj) + /iQL/3i(n]72 • 0* • J.Pj • u;*)] )^ 



It should be clear that many factors (i.e. basis vectors w"^, , v and operators P+) can be 
taken out of the correlators as they do not depend on the integration variables («-, /?+, 

(p) but are simply functions of momenta; we do not explicitly do this here for brevity. 

A few comments are now in order. The measure current contribution exactly cancels the 
similar contribution of the light sector, the first term in (3.17), and, as stated many times, all 
the "action" is in the other terms. Consider first the terms independent on k. The correlators 
that enter Ufiu (4.12) scale as 1/y and l/y^, respectively, so their contributions to 11"^, are 
y- independent. As for the /i-dependence of the large-y limit of n~^, we only know for sure 
that h should drop out of the divergence of the polarization operator (it is difficult to see how 
this can be proved starting from the explicit expression of the mirror polarization operator, 
but it must be true, as a consequence of the gauge invariance of the full "light" + "mirror" 
partition function and the local smoothness of the split). Finally, we expect that the terms 
proportional to k, 11^^ and H^u vanish in the k — > limit: it is not clear how a 1/k singularity 
in the scalar Green's functions could come about if there is indeed a disordered phase in this 
limit (as per the results of [38] ) ; we numerically test the ^ behavior of II^j^ + YlJ^u in 
Section 4.2.5 and confirm this expectation. 

In our numerical simulation, we compute II^j^ of eqn. (4.11) and look for non-local terms 
in its real part. Before stating the results of the simulations in the next Section, we note 
that in order to compute the scalar and fermion correlators in (4.11), one needs to choose 
a representation of 7-matrices, a basis of 75-eigenvectors, work out the perturbative expan- 
sion of the Neuberger-Dirac operator to second order in the gauge field, and finally, express 
the correlators (4.12-4.14) in terms of these quantities. We give our notation and conven- 
tions in Appendices B and C in sufficient detail to complete this straightforward but tedious 
calculation. 



4.2 Monte Carlo results for the mirror polarization operator at strong coupling 

In the continuum, the contribution to the Fourier transform of the real part of the polarization 
operator due to massless particles is: 

%,{k) = 2C '^M^^'^- k^k, ^ 2CgS scalar ^ ^Cch.ferm. - " -0-159, (4.15) 

where the constant C depends on the number of massless degrees of freedom and the normal- 
ization is chosen for later convenience. The values of C given are for a single Green-Schwarz 
scalar and a single chiral fermion, inferred by comparison with (3.24) and (3.23). The polar- 



- 26 - 



41,1] Tll-21 



0.05 - 



-0.05 - 



-0.10 - 




Figure 1: The real parts of Hn (left plot) and (right plot), the Fourier components of the real 
part of the basis-vector independent part of polarization operator of a free chiral GW ferniion on an 
8^ lattice, as a function of momentum. The /c-dependence is plotted along three lines in momentum 
space — the circles and horizontal squares correspond to lines approaching the origin at O" {k2 = 0) and 
45" (fci ~ k2 ^ k) wrt x-axis, respectively, while the tilted squares are taken on a line with fci = (at 
90"). Clearly, the discontinuity at small k (i.e., at k — 7 and fc = 1) is as predicted by the continuum 
expression — for example, the value of 2C read off nii(90°) agrees well with the continuum expression. 
We note that in these plots (generated by Mathematica) , momenta are labeled by 1...8 instead of 
0...7 as in the plots presenting the Monte-Carlo results; an identification = 8 should be made for 
comparison. 



ization operator (4.15) has a directional singularity as A; — > 0: 

nii(0)L_,o = C{l-cos2cP) 

n2i(0)L_o = ' (4.16) 

where (p is the angle of approach to the origin measured from the positive-fei axis. Thus, from 
(4.16) we expect that on the lattice, if there are massless particles, the following relations will 
hold as /c — > 0: 

-n2i(45°) = c , 

2C , (4.17) 

n2i(o°) = n2i(90°) = . 

To test the relations (4.17) on the lattice, we begin by first computing the real part of the 
polarization operator of free massless GW fermions, which is simply equal to 1/2 the Dirac 
GW fermion polarization operator and the analytical expression for which is given in (A. 5) of 
Appendix A. On Fig. 1, we plot the real part of the Fourier transform of the real part^^ of TE^^ 
for a free chiral GW fermion (A. 5); the Fourier transform is precisely defined in D.2. In order 
to make the calculation of (A. 5) well-defined, we imposed antiperiodic boundary conditions 

^'^This cumbersome expression is, unfortunately, unavoidable: real and imaginary parts of 11^,^ are defined 
in x-space, but then the finite lattice Fourier transform (D.2) of Rell^^ has both a real and imaginary part 
and becomes real only in the continuum limit. 



nii(45°) = 
nii(9oO) = 
nii(o°) = 
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on the fermions, which results in a well-defined value at fc = 0, but does not remove the 
directional singularity as A; — > 0. 

It is clear from Fig. 1 that the polarization operator approaches different limits as one 
approaches the origin at different angles. The discontinuity at small momenta is clearly 
visible already on the rather small lattice used. The results agree well with the expectation 
of eqn. (4.17). The value of C inferred from A: = 1 on an 8^ lattice does approximately match 
the one of the continuum result for a chiral fermion, given in (4.15); a rather precise numerical 
agreement can be seen on somewhat larger lattices (e.g., virtually identical on 32^), however, 
we only plot the 8^ result since this is the lattice size of our numerical simulation. 

If one performs a similar calculation for a massive fermion, the /c — > discontinuity 
disappears, in accord with the expectation that a massive-particle loop contributes 11^,^ ~ 
^ ^('^/.ti/A;^ k^kjy). 

In what follows, we will compare the results for Tin and 1121 for the free chiral GW 
fermion with the results of the Monte Carlo simulation of the same components of the mirror 
polarization operator. While in the text we refer to 11 and 12 components of the polarization 
operator, in the Figures and the captions, the replacement IIn IIoo, II21 IIio, will be 
used. 

4.2.1 A few v^ords about the simulation 

All our simulations are on an 8 x 8 lattice. The reason we only consider such a small lattice is 
that the computation of 11^^ is rather demanding, because of the large number of momentum 
sums that occur in the correlators, most notably in the fermion-fermion current-current cor- 
relator of eqn. (4.11). Nevertheless, when combined with the analytic results of the previous 
Sections and the Appendices, even this rather small lattice is sufficient for a qualitative study 
of the mirror spectrum and in particular, an identification of the massless mirror states. 

The code used to generate configurations was developed by J. Giedt for the study of [38] 
and uses the cluster algorithm to generate AY-model configurations with a reweighting of 
the fermion determinant. The fermion measure is taken into account through determinant 
reweighting: 



Here, O is any observable, and (• • ■)r^ denotes an expectation value with respect to the measure 
of the XY model (cf. eqn. (2.2) with U = l), 



We monitor the reliability of this method in several ways. We measure the autocorrelation 
time for reweighted quantities (OdetM)^, as well as (detM)^, to be certain that the config- 
urations remain independent with respect to the new measure. We perform a jackknife error 
analysis of the averages {O) that are obtained, gathering sufficient data to keep errors small. 



{O det M)n 
(det M)rj 



(4.18) 




(4.19) 
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Figure 2: The real parts of IIoo and IIio of the mirror for k = 0.1, = 2, as a function of momen- 
tum approaching the origin in different directions: symmetric phase, massless fermion (compare with 
Fig. 1.). 
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Figure 3: The real parts of Hon and IIio of the mkror for n = 0.1, h = 5: symmetric phase, massless 
fermion (compare with Fig. 1.). 



An additional package was developed by one of us (Y.S.) so that the original code was 
adapted to work on the Sunnyvale computer cluster at CITA, where the momentum sums 
over different configurations were performed with hundreds of processors in parallel. In order 
to reduce the statistical errors, 16000 independent configurations were generated and used in 
the calculation of the averages. 

4.2.2 Strong-coupling symmetric phase {h> 1,k <^ Hc)'- massless chiral fermion 

We now compare Fig. 1 to the numerical simulation of 11^^ of the mirror theory, an expression 
for which in terms of mirror correlation functions is given in (4.11). On Figs. 2 and 3 we 
show the same components of the polarization operator as on Fig. 1 for k = 0.1 and h = 2,5, 
respectively. On Figs. 4 and 5 we show the results for k, = 0.5 and /i = 2,5, respectively. In 
all cases, measurements of the susceptibilities [38] indicated that the theory is in the strong 
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Figure 4: The real parts of IIoo and IIio of the mirror for k — 0.5, h — 2: symmetric phase, massless 
fermion (compare with Fig. 1.). 




Figure 5: The real parts of IIoo and IIio of the mirror for k = 0.5, h = 5: symmetric phase, massless 
fermion (compare with Fig. 1.). 



coupling symmetric phase. 

Thus, we observe that for k < 1 and h > 1, i.e. in the strong-coupling symmetric phase 
found in [38] , the polarization operator of the mirror theory is qualitatively and quantitatively 
close to the one in the free GW-theory of Fig. 1. The discontinuity at small- A; is as expected 
from the continuum formula for the real part of a chiral fermion polarization operator and the 
actual numerical values of Hell^u are also close to the ones for the free fermion (for example, 
the values of 2C inferred from nii(90''), see (4.17), are closer to the continuum value for 
K = 0.1, deeper in the symmetric phase, than for k, = 0.5). The real part shows slight 
variations with h, as the plots for = 2 and h = 5 show, as well as with k, (when k < 1), as 
the plots for k = 0.5 and k = 0.1 show. In each case, the ratio of nii(90'^)/nii(45*') ~ 2.5, 
are in agreement with the ratio for free GW fermions of Fig. 1 and somewhat larger than the 
continuum ratio of 2. 

In every case, the divergence of the imaginary part, as well as eqns. (3.27), and (3.29), 
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Figure 6: The real parts of IIoo and IIio of the mirror for k = 5, /i = 2: broken phase, massless 
scalar. 

were numerically checked to be obeyed by U'^,^ up to the errors of the simulation. We conclude 
that in this phase the mirror theory satisfies 't Hooft anomaly matching in the massless chiral 
fermion mode. 

4.2.3 "Broken" phase (h > 1, Kc): Green-Schwarz scalar 

The results for the mirror polarization operator in the broken phase are shown on Fig. 6, 
for K = 5,h = 2. Taking k ^ 1 takes the unitary Higgs field into an algebraically ordered 
( "broken" ) phase. The studies of [38] showed that the results of simulations in this phase agree 
very well with perturbation theory (a ^ expansion) . This continues to be the case here — note 
that the values of at small momenta scale with k, in contrast to the symmetric k < 1 
phases, see eqn. (3.23). From nii(90'^) we infer 2C ~ —5, in good agreement with the value 
for K = 5 of eqn. (4.15). Thus, we conclude that the discontinuity of 11^^ on Fig. 6 in the 
"broken" phase is due to a massless scalar. 

4.2.4 Strong-coupling symmetric phase (/i — > 0, k <^ Kc): massless chiral fermion 
and massless vectorlike pair 

At /i = 0, in Section 5, we give analytical arguments that the massless spectrum should 
contain three charged massless chiral fermions (two of them forming a vectorlike pair and 
thus not contributing to the anomaly). This is borne out by the simulations done in the 
h ^ limit. For example, the value of |2C| inferred from the plot on Fig. 7 of nii(90'^) 
is roughly .39, while that inferred from 1112(45'^) it is .44 (closer to the continuum value of 
3/(27r) ~ .48). A calculation (which for brevity we do not show) of the polarization operator 
of the mirror partition function det+(2 — D), see eqn. (5.6) in Section 5, also shows good 
agreement with the h ^ Monte Carlo data shown on Fig. 7. 

4.2.5 Scaling of the different contributions to II^i, as k — > 

Finally, on Fig. 8 we show the K-scaling of the K-dependent terms in 11^^ of the mirror 
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Figure 7: The real parts of IIoo and IIio of the mirror for k = 0.5, h = 0.01: symmetric phase, three 
massless chiral fermions, see Section 5; the dip at fc = 4 is due to an "Umklapp" process (an analytic 
calculation using the mirror partition function det+(2 — D) gives a very similar result). 



theory (i.e. deep into the symmetric phase). More precisely, we show the 0-90 degree split 
of the sum of the IIqq (4.13) and IIqq (4.14) at small momenta as a function of k. It is 
clear that the contribution of the K-dependent terms scales to zero with k, showing that in 
the symmetric phase contributions to the angular discontinuity of polarization operator come 

2 

from the Yukawa terms, IV'^ of eqn. (4.11). 



5. Attempt at (some) analytic understanding 

In this Section, we will attempt to get a somewhat better analytic understanding of the nu- 
merical results of the previous Sections. To this end, we will develop a different representation 
of the 1-0 model partition function. We begin by noting that the 1-0 model action (2.1) can be 
written in unprojected components (the superscripts P^^^^ in the hatted projectors indicate 
whether the charge-1 or charge-0 Dirac operator has to be used): 

Skinetic = - (ipDilp) - ixDoX) 
SYukan,a = y { (V^Pj^^P+x) + (xP-<t^P~^) + h[{^'^ {P^f (t^l2P+x) " [xP-l2^* {Plf^^)] 

where the measure now is the usual vectorlike theory measure in terms of -0, V', X) X- 

Z = /*#<iX«#.-^— (5.1) 
and S'k is defined in (2.2). We now perform the field redefinition: 



Splitting of P\qq^^^^' +P\qq'' ^^^^ near zero momentum, h=2 
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Figure 8: The small-K scaling of (Jlli + Il\i)\^o of the K-dependent terms' contribution to the 
discontinuity of the mirror polarization operator at small momentum. 



The motivation for this redefinition can be traced back to the GW relation, which implies 
that: 

^ 75 = 75 ^ , (5.3) 

in other words, a 75 action on the original tp is transformed into the action of 75 on the new 
fields in (5.2). The price to pay for having an action of the lattice chiral symmetries generated 
by 75, exactly as in the continuum, is the nonlocality of the redefinition (5.2) and thus of the 
resulting action (see (5.4) below). 

Since we work perturbatively in the gauge field and at finite volume, we will imagine 
throughout this Section that the singularity at = 2 of (5.2) is avoided by turning on 
background Wilson lines. The effect of the field redefinition on the partition function is: 

S' —S' 

u^^^ + lf2 Yukawa 



Z = y di;d'iljdxdxd(t> det(2 - Di) det(2 - Dq) e~^k^net 
S\.^etrc = - ^) - (x x) , (5.4) 

+ yh (V'^<^72^+X) - \ (xP-i2<P*i^^) - \ (x^-^^72</'*(^^)^^'+V^^' 
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Obtaining the transformed Yukawa couplings in "Sy^^^^^ requires repeated use of the GW 
relation and the equivalent relation (5.3). Needless to say, similar redefinitions hold in four 
dimensions and can be straightforwardly performed if necessary. 

Several comments, concerning the representation (5.4) of the mirror partition function, 
are now in order. We hope that these comments are useful to clarify the relation between 
different formulations of exact lattice chirality: 

• The singularity of the action at the position of the doublers, D = 2, ensures that they 
have infinite action and do not propagate at the classical level, as first proposed by 
Rebbi [52]. The problem with [52], pointed out in [53, 54], of the doublers contributing 
as ghosts to the photon polarization operator at the quantum level is solved by the 
determinant prefactors, which exactly cancel the would-be ghost /doubler contributions. 

• Another comment concerns the relation of (5.4), with y = /i = 0, to domain wall 
fermions. In the 2d case, these propagate on a finite interval in three dimensions. The 
generating functional of Green's functions of the boundary chiral modes — in the case 
at hand, one charged and one neutral — can be obtained by integrating out the bulk 
fermions. This is technically possible since the fermion action is bilinear and there is 
no gauge field propagation in the extra dimension (or any other non-uniformity except 
at the boundaries). When taking the chirally-symmetric limit of an infinite number of 
sites in the extra dimension, a massive bulk Pauli-Villars field, antiperiodic in the extra 
dimension, has to be included in order to obtain a finite determinant. The generating 
functional of Green's function for the boundary chiral modes can be represented as a 
partition function with source terms. The result is exactly (5.4), with y = h = 0, and 
with source terms for ip and x included. The determinant prefactors in (5.4) arise as a 
combination of the determinants of the bulk fermions and the Pauli-Villars fields. The 
derivation of these results can be extracted from the work of ref. [56] (a less rigorous but 
nonetheless enlightening continuum discussion of the properties of the kinetic operator 
for the boundary modes is given in [57]). 

• For vanishing Majorana coupling, h = 0, the Yukawa interaction in (5.4) is equivalent 
to the chirally invariant Yukawa coupling of ref. [44], see also [58]. To see this, use 
(V'j^V') = ^{'^Dtp) + ^{i/jD^^i^) to replace the ip kinetic term in the action S^.^^. 
Then note that -\{^D^ij) = (^(2 - D)C) + ^(^-D?/') + -^^DS), where it is under- 
stood that the new charged field is integrated out from the action using its equation 
of motion, while the ^-determinant exactly cancels the one in (5.4). Then, shift the 
integration variable ip ^ ip + Next, perform exactly the same operations on Xi 
introducing a neutral field rj, to finally obtain the action in the original form of Liischer, 
suitably adapted to the 2d case and to our normalization: 

(5.5) 
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where the measure is the trivial one over tp^Xi'H^i- The Yukawa interaction of (2.1) is 
thus equivalent to that of [44]. We note that gauge and chiral invariant Majorana cou- 
plings were not considered in [44] and to the best of our knowledge were first constructed 
in [38]. 

We can now use the representation (5.4) to split the partition function exactly as we did 
in the basis of GW fermions. The representation of the split partition function that we give 
below is, in fact, equivalent to that in (2.11). The splitting of the ijj-x partition function (5.4) 
into a "light" and "mirror" part is now done trivially using the 75-eigenvectors, which have 
no gauge- field dependence. We obtain, denoting now by V^iiXi the "normal" 75-chirality 
components of the 2-component V'jX) and using the fact that Syukawa o^^y depends on the 
mirror components of V') X'- 

Z = Z+ X Z- X J , 

Z+ = det_(2 - L>o)det+(2 - Di) J di}+di}j^dx^dx^e~^+'^'^+)~^-'^^~) , (5.6) 

Z_ = det+(2-L>o)det_(2-Di) j d'0_#_dx+dx+#e~('^"^'^")~(^+^^+)"^^"'=''-""^'' . 

Splitting the determinant prefactor into "light" and "mirror," as indicated in (5.6), requires 
using the gauge- field dependent eigenvectors of 75 : 

det+(2-Z)) = det+(l + 7575) = det||(n|(2-D)t,-)|| =det||2(^xjt,)ll, (5-7) 
det_(2 -D) = det_(l + 7575) = det||(u;|(2 - D)vj)\\ = det||2(?i;Jt;j)||, 

where the appropriate 75-eigenvectors are to be used for or . The gauge variation 
of the mirror partition function Z_ has now two contributions: one from the variation of 
det_(2 — D^) and one from the variation of the path integral over ■ij;-,x+,4'- The splitting 
theorem, applied to the chiral partition function defined by the latter, can be easily seen to 
imply that the gauge variation Z_ with the determinant factors left out vanishes. Thus the 
entire gauge variation of Z^ comes from the determinants: 

<51nZ_ = <5^1ndet||(^i;J(2 - D^)v.i)\\ = J^(5a;U;]u;,) + tr P+^-i^<^^(2 - D^) 

j 

= J^(5^^«;,) + iTiu (P| - P+) = + '-Ttuj%, (5.8) 

j j 

where, as usual, the measure current is cancelled by the variation of the Jacobian and light 
partition function. The gauge variation of Z_ is, naturally, the same as in (2.15) (the gauge 
variation of Z+ can be obtained similarly to (5.8) and be seen to combine, together with the 
Jacobian to cancel that of Z^). 

The determinant prefactor in Z^ contributes both to the real and imaginary parts of 
the mirror polarization operator. It is clear from (5.8) above that the imaginary part of the 
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polarization operator due to the determinant is exactly as required by anomaly cancellation. 
Assuming unitarity at long distances, one would argue that Rell^jy should receive a contribu- 
tion from at least one massless charged chiral fermion (plus, possibly, a number of massless 
states in anomaly-free representations). The goal of our simulation was to precisely find out 
the real part of the mirror polarization operator. Let us now compare the numerical findings 
with what can be inferred from the representation of the mirror partition function (5.6, 5.7). 

1. The real part of the mirror theory 11"^ receives two contributions at y = oo. The first is 
due to the determinant prefactor. The contribution to Ren"^, of det_(2 — D^) is easily 
seen to be exactly that of three massless propagating charged chiral fermions — the chiral 
components of the three 2d doubler modes. This can be easily done using the formulae 
from the Appendix to calculate the determinant contribution to n^jy, then numerically 
plotting the result and comparing to the long-distance contribution of a single charged 
chiral fermion. The determinant prefactor contribution to Ren~j^ is, obviously, the same 
for any value y,h. 

2. The only other contribution to Rell"^ at y = oo arises from the nonlocal coupling in 
the Majorana Yukawa term in Sy^^j^awa i^-^)- While it appears difficult to calculate 
analytically this contribution in the disordered- phase, our numerical results for the 
mirror polarization operator Rell^jy show that there is one massless propagating chiral 
fermion. Thus two of the three massless modes contributed by the prefactor are can- 
celled by the nonlocal contribution to Rell^i, from the Majorana term in Syukawa] it 
appears that this cancellation is exact for all values of /i > 1. 

3. For vanishing Majorana coupling, /i = 0, on the other hand, both numerical simulations 
(with the code of [38] used in this paper, we can only approach the /i ^ limit) 
and analytic arguments using the representation (5.6), to be discussed in more detail 
elsewhere, show that the real part of the mirror polarization operator at small momenta 
and infinite y is that of three charged massless chiral fermions. This is not too surprising 
as the Majorana coupling goes away in the limit and at infinite y the only dependence on 
the gauge field is in the determinant prefactor in Z_ of (5.7). This indicates the crucial 
role of the Majorana-type couplings (recall that they were motivated by the need to 
break all mirror global symmetries) in facilitating the decoupling of the maximal possible 
number (allowed by anomaly matching) of charged mirror degrees of freedom. 

4. Another lesson we learned is that probing the fermion spectrum with local fermion op- 
erators (including charged local fermion-scalar composites, as in [38]) can miss massless 
degrees of freedom. The massless charged mirror fermions were not seen in that study, 
perhaps because they are not expressed in an obvious local way through the original 
variables. The long-distance gauge-boson polarization operator of the mirror theory is 
a universal probe of the charged mirror spectrum and should be the first quantity, along 
with susceptibilities probing chiral symmetry breaking, computed in any future studies 
of anomaly- free models. 
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A. More fun with chiral polarization operators 

Consider first the mirror polarization operator (4.7) at y = 0, which can be written as 
= 5yj^ + n^. The part of the polarization operator is independent of the basis 
vectors and can be expressed (easily verified by looking at (4.7)) as: 

n'± = s^-W^ , (A.i) 

with: 

n'± = Ti'P±D~'S^D = ^TiD-'6^D ± ^Trj^D^'S^D . (A.2) 

We will denote the two terms in (A.2) as the vector (superscript V) and axial (superscript 5) 
parts of n^: 

Ul = Ti^^D-^d^D . (A.3) 

We note that Iljf is real, since = 75D75. On the other hand, we easily obtain (using the 
GW relation on the way) that: 

Ren^ = ^TY75<5^Z) 
ilmn^ = iTY[(^^A75]^-' . (A.4) 

Thus, we can express 

n^t (A.I) as follows: 
n'± = \ S,TtD-^6^D ± J Tt-,^S^5,D ± ^ 6,Tt [6^D, 75] (A.5) 
A few comments: 

1. The first two terms of 11^ are real and symmetric in /d, v, while the last term is purely 
imaginary and has both a symmetric and an antisymmetric part. 

2. The first term, ^5i,'TtD~^6^D, in (A.5) is simply 1/2 the vector theory polarization 
operator. It is manifestly real and has a nonlocality (a factor of D~^) due to the fact 
that a massless particle is propagating in the loop. It is also manifestly transverse w.r.t. 
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both indices. In the continuum hmit it reduces to the transverse Lorentz invariant 
contribution (the first term in (3.24)) with a coefficient equal to one-half that of the 
massless Schwinger model. 

3. The second term, ±^TT:j56fj,6uD, is manifestly symmetric and is local (with the usual 
exponential tail). It is transverse w.r.t. both and i'. This term should vanish in the 
continuum limit, as there is no symmetric, local, transverse, parity-odd expression one 
can write in the continuum. In fact, numerical evaluation of the trace on finite lattices 
indicates that this term vanishes identically at A = 0; we have no analytic proof, but 
feel that one should exist. 

4. The last term, zt^Jj^Tr [6^D, 7,5] D~^, is purely imaginary, and has no manifest symmetry 
in /i, u. We know from general arguments that its antisymmetric part should equal 

(see also (A. 17) below) and thus be local. Its symmetric part, on the other hand, is 
nonlocal. The divergences of the symmetric and antisymmetric parts are equal, each 
giving rise to one-half the anomalous divergence 

ofn'^t- Once again, this is easy to also 

explicitly check from (A. 5): 

V*Xt{x, y) = ±\vl <5,TY [5,D, 75] D^^ = U^, x) . (A.6) 

Two questions are to be answered in order to generalize these results to the most general 
chiral theories (e.g. to our y 7^ mirror). Our object of interest 11"^^' is again defined by: 

n;. = 5,i- + n;;. (A.7) 

First, is it true that Il~J is always a total derivative? Namely, does there always exist a n~ ' 
such that: 

Second, if such a n~ ' does exist, is it manifestly gauge invariant? 

The answer to the ffist question is immediately yes if one just looks at the definition 
of n^j^' given by equation (A.7) and recalls that both II^j^ = 5^(5^ log Z_ and S^j^ are total 
derivatives. In fact, we know precisely what 11" ' is. Given any chiral action, in the notation 
of [41]: 

S[X,Y\0]^eMS[X,Y\0]), 

and the partition function: 

Z = j UdcdcS[ciUi,Civl,0] , (A.8) 
with Ui and Vi some appropriate eigenvectors, we have proved the "splitting-theorem": 

7(1) 

6^logZ = j^ + ^, (A.9) 
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where: 

Z«= Udcdc6o,^,S[cm,c,vj,0]. (A.IO) 

Here So,fj,S represents the variation of S with respect to Afj, with the vectors Ui and Vi fixed 
(namely, the variation of S due to the variations of the operators only). As explained in 
Section 4.1.1, to make sure 5o,ij.S[ciUi,Civl, O] is again chiral, one should insert additional P's 
if needed, and the value of zj^^ remains intact in this procedure. We will implicitly assume 
this rule here. Comparing equation (A. 9) with equation (A. 7) implies: 

n;; = 5j^j. (A.ll) 

Now, the answer to the second question is also yes provided that S is gauge invariant. 
We assume that the action S'[X, y^^,0] is invariant under the gauge transformation 

5X = iujX, 6Y = iujY, and 60 = i[uj , O]. (A.12) 

It is easily verified that Sq.^iS is also gauge invariant, where all the operators, including 5^iO 
and the additional P's inserted, transform as adjoints as well. Therefore both Z and zj^^ are 
partition functions defined by a "chiral" and gauge invariant action. As given by equation 
(4.55) in [41], their gauge transformations are completely determined independent to the 
actual expressions of the actions. Indeed, 

6^ log Z = 6^ log = - iTrcuiP - P) , (A.13) 

where ji^ is the measure current corresponding to an infinitesimal gauge transformation with 
parameter lo. Thus, we find: 

z(i) 

6^ log -^ = S.o log - 6^ log Z = 0, (A.14) 

which says that log is gauge invariant, and so is ^y~- 

Finally, we derive the general properties of chiral polarization operators listed in Section 
3.4. Gauge invariance of n~ ' implies immediately, that in a general classically gauge invariant 
chiral theory: 

V:U-J = 0, (A.15) 

following the steps that led from (3.11) to (3.14). Furthermore, in Section 3 we found that 
the divergence of with respect to the first index fi is independent of the action: 

V;n;; = '-6.tTjl. (A.16) 

Another action-independent result is that the anti-symmetric part of Il~J is also known 
explicitly — since is manifestly symmetric, the antisymmetric parts of 11"^^' and S^jj^ 
should cancel: 

n^; ^ = ^i^'^J - n;;) = -:f^., (a.i7) 
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v;;n;;^ = ^v*n;; (A.20) 



where, 

T^, = <5^jr - = -Tr(P-~ [ , ] ) (A.18) 

is manifestly local. Its divergence, after a few steps of algebra, turns out to be ^S^tT^xx ^ 
well. Applying (3.13) to .F^j/ and using the gauge transformation d^jP- = itoP- — iP-to, we 
find: 18 

r 

= i —^t[5uP-P-{ujP- - P-lo) - P.6^P.{ujP- - P_a;)] 

6 . ^ (^-19) 

Combining these results leads to: 

2 

In other words, the anti-symmetric and symmetric parts of 11"^' each contribute half of the 
anomalous divergence. The results (A. 15), (A. 16), (A.20) just derived hold for general chiral 
partition functions, independent of the action, and in particular for our y mirror theory. 
Verifying that these exact properties hold is an important check on any numerical simulation. 

B. The Neuberger-Dirac operator and its perturbative expansion 

The Wilson operator in our convention is given by: 

A* 

+ ^ ^ {^5rnn - 5rn+fi,nU^{m) - 5rn,n+ilUl[n)^ , (B.l) 

where m,n label d-dimensional lattice sites, U^{m) = Q^^i^i"^) ^ /i is a unit vector in the //-th 
direction on the lattice, and 5m,n is a d-dimensional Kronecker symbol. The Wilson operator 
(B.l) is 75 hermitean {X'^)mn = l^XmnH^- We define 75 as follows: 

in terms of which our convention for the GW operator is: 

D = l- 7575 = Dq + Di+D2 + ... , (B.3) 

where -Do, 1,2 are the terms in the expansion of D around the trivial gauge background. We 
take M = 1, r = —1 from now on. We define our finite-volume Fourier transforms as follows, 
using cvn = e n : 

X{q,p) = -^Y.^'n'"''''"'^^- ■ (B.4) 



^We used identities P^SP-P-. = and 5P-P- + PSP- = 5P-, both from PI 
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We also define the functions: 



Tip Trp 
s{p) = sm — , c{p) = cos — . (B.5) 

By (To, we denote tlie 2 2 -dimensional unit matrix; in d=2, we use 71 = cr^, 72 = 0"^,75 = cr^, 
where a* are the Pauli matrices. Similar to (B.3), we also have an expansion of X in powers 
of Af,: 

X{p,q) = Xo{p,q) + X^{p,q) + X2{p,q) + . . . , (B.6) 
where the trivial-background term is given by: 

Xo{p,q) = Sp,qXo{p), where xo{p) = I 1 - ^2s^(p^) | ctq - ^ij^s{2p^) . (B.7) 

\ M / M 

The linear and quadratic terms of (B.6) are: 

Xi{p, (l) = -Jpl^ ^t^{P -q)^n'' [«7/^c(Pm + 9m) + cros(PM + , (B.^ 

^2(p,g) = -^E 



where in Xi we use the Fourier transform of A^, defined by: 

A,{p) = Y,u;-/-'^A,{7n). (B.9) 

m 

In terms of Xq, Xi, X2 the various terms in the expansion of the GW operator (B.3) are 
described below. We begin with the free GW operator, specializing to two dimensions in the 
second line below: 

Do{p,q) = 1 , ^0 = Sp,q do{p) , 

where, consistent with the notation of [38] , we have defined the function of momenta: 

1 - 2s(pi)2 - 2s(p2)' 



a{p) = 1 



w{p) 

,s(2pi) 



m ^ '-^ . (BAD 

w[p) 

d{p) 



w{p) ' 

w{p) = v/l + 8s(pi)2s(p2)2 . (B.12) 
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The linear term in (B.3) is given by: 

Di{p,q) = (-X,{p,q) + ^ Xl{p,q) ^) , (B.13) 

w[p) + w[q) \ w{p) w[q) J 

where xl{p, q) = 75X1(^5, (7)75, and w, xq, Xi are as defined above. The second-order term in 
(B.3) is then found to be: 



w{p) + w{q) \ ' w{p) ' w{q) 

^ 1 
s^^x [^(^') + + w{k)][w{q) + w{k)] ^ (B-14) 

{fcM = l| 



+ 



w[p)w[q)w[k) 

xoip) X\{p,k) Xi{k,q)+Xi{p,k) X\{k,q) xo{q) + Xi{p,k) xj(fc) Xi(A:,(?)} . 
Both the hnear and quadratic terms in the expansion are consistent with the ones found in 
[55], the precise map being £)^^^p^p^^ = _ ^D^^2^^^^^ 



, usmg our expressions 



for xo,Xi,X2. 
C. The chiral eigenvector basis 

To compute the mirror polarization operator, we need explicit expressions for the expansion 
(2.9) in terms of 75 eigenvectors for ^ = 0. These were worked out in the Appendix of [38] 
and we give the A = expansions (2.9) here: 

X+ix) = ^l3lvk{x) , x+i^) = ^f_,ul{x) (C.l) 

k k 
k k 

where the sum is over momenta, k\,k2 = 1 . . . N, which label the independent eigenvectors 
for a trivial gauge background. 

The 75 (75) eigenvectors t, v {w, u) are orthonormal (we note that the phase factor in 
tk{x) is a matter of pure convenience and is included for agreement with [38]) and are: 



Vkix) 

4(^) = ^^^'''(0 1)^fc, (C.2) 



tk{x) 
wl{x) 



1 




N 




1 




N 


^^'■^ ( 


1 




N 


w^^' e" 


1 


-k-x I 


N 
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Here we use the unitary matrices C/^ and 14: 

2 -Xle^^" 2-Xl J ' 2 y Afc (2 - Afc)e-^^'= ^ ' ^ ' ^ 

defined in terms of the positive imaginary part eigenvalues of the GW operator (the functions 
a, b, d are as in (B.ll)): 

\k = a{k) + i ^h{kY + d{kY , (C.4) 

and the phase factor: 

e^.. ^ / tUSw' (^-^^) ^ iN,NUN/2,N/2UN/2,N),iN,N/2). 
\ 1, (A;i,A:2) = (Af,Af),(iV/2,iV/2),(iV/2,iV),(iV,iV/2) . 

Using these definitions, it is straightforward to verify that the w, u and t, v bases obey 
Sa:(^I(^) ■ ^p(^)) ~ ^kp, J2x('^k(^) ' ^p(^)) ~ 0, etc; proving this only requires (in the 
w,u case) use of the GW relation + A|. = |Afcp. 



D. Perturbative derivation of the anomaly 

Here, we derive eqn. (3.18). Now, we know that — see eqn. (3.17), as well as (A. 5, A. 6): 

i 6tic %{x,x)[A] 



(D.l) 



A=0 



For the purpose of comparing with finite-volume simulations of the mirror sector, we would 
like to have a finite volume expression for the r.h.s. of (D.l). 

To begin, we first define the Fourier transform of any free (i.e. translationally invariant) 
polarization operator H^^-uy similar to the usual continuum definition: 



(D.2) 



We then also define the Fourier transform of its divergence: 



(D.3) 



Finally, we similarly define the Fourier transform of the r.h.s. of (D.l): 

.k.y Jtr 75(0,0) [A] 



tu{k) 



E 



N 



A=0 



(D.4) 



allowing us to recast (D.l) in the form d^{k) = ^ty{k). 
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Now we can compute the function tu{k) explicitly, using (B.3) and (B.13) for Di, and 

find: 



N 



" ^ N\^^ w{p)w{p - k)[w{p) + w{p - k)] ^ 

{^2^, - k,) e^^Uv)UP -k)- c{2p, - k,) e^^ {i^{p) - ^p - k)]] , (D.6) 
where e^^ = 1 and we have defined: 

f{p) = 1 - 2s{p,f - 2s{p2? (D.7) 

It is easy to check that in the small-A: limit, eqn. (D.5) reduces to: 

, 27Tk^J_ A c(2pi) + c(2p2) - c(2pi)c(2p2) 

tu[k)^eu^ iV 2iV2 2^ w{p)^ 

{Pa=l} 



+ o{e) (D. 



where A;^°"*- = is the continuum momentum; we note that this establishes the coefficient 
in (3.18). 
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